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We generalize the nonlinear sigma model treatment of quantum spin chains to cases including 
' ferromagnetic bonds. When these bonds are strong enough, the classical ground state is no longer 

, the standard Neel order and we present an extension of the known formalism to deal with this 

■ situation. We study the alternating ferromagnetic-antiferromagnetic spin chain introduced by Hida. 

The smooth crossover between decoupled dimers and the Haldane phase is semi-quantitatively re- 
> | , produced. We study also a spin ladder with diagonal exchange couplings that interpolates between 

O i 1 the gapped phase of the two-leg spin ladder and the Haldane phase. Here again we show that there 

| is a good agreement between DMRG data and our analytical results. 

m ■ PACS: 75.10.Jm, 75.20.Hr 



D 1 I. INTRODUCTION 

h. 

c/2 , One-dimensional quantum spin systems exhibit remarkable physical properties. One of the most interesting case is 
the s=l -antiferromagnetic Heisenberg spin chain. Contrary to the s=l/2 case, this system has a gap as predicted by 
HaldaneEJ and a finite spin correlation length. It is an example of a system which is disordered at zero temperature 
due to quantum fluctuations. The original conjecture has been checked experimentallycl, numericallyu as well as 
analyticallyDtl. Although the,ground state is disordered in the sense that spin correlations decay exponentially, there 
is a hidden topological orderQ that is revealed in the bulk of the chain only by nonlocal observables or by ground 
state degeneracy in an open geometry. This hidden order is most clearly seen in the VBS wavefunction which is an 
O [ approximate ground state of the spin-1 chainO. To construct this wavefunction one has first to write each spin s=l 
as a triplet of two fictitious spins s=l/2. Then one couples nearest-neighbor spins s=l/2 into singlets. This leads to 
| a function which is obviously singlet and translation invariant. It is an excellent approximation of the true ground 
. stateB. The perfect crystalline pattern of singlets is the hidden order. K. Hida has given an appealing picture of 
0^ ' the Haldane gap and the VBS ground state by considering an alternating ferro-antiferromagnetic s=l/2 chainO. The 
Hamiltonian is given by (see figure 1) : 



S 

i 



O ■ 

on : 



c3 



X 



H. — Jaf J]] &2n ■ S2n+1 + Jf J]] S2„+l • S2n+2- (1-1) 



Here are s=l/2 spin operators, Jaf is positive and Jf is negative. In what follows, we set Jaf = 1 and Jf = — 7- 
The family of systems defined by Eq.([hl]) has simple limiting cases. For 7 = we have a set of decoupled pairs of 
spins that have a trivial ground state : all pairs are locked in singlets. When 7 — *■ 00, the s=l/2 are coupled by pairs 
into s=l states and we get a chain of spins s=l. Hida has studied numerically the gap of the system as a function of 7. 
He has used Lanczos diagonalization techniques to evaluate the gap and he showed that there is no phase transition 
as a function of 7. As a consequence, the Haldane gap of the limit 7 — > 00 is continuously connected to the trivial 
gap of the decoupled limit 7 = 0. The excitation spectrum also evolves smoothly. So the alternating chain offers a 
simple physical picture of both the Haldane gap and the hidden topological order. 

There is another approach to quantum spin chains which is the continuum field theory known .as the nonlinear 
a model (NLtrM). Introduced originally by HaldaneEJ, this field theory includes a topological tern£3 9 — 2ns for a 



spin-s chain. While this term can be discarded when s is integer, it is responsible for masslessness when 9 — tt (mod 
27t), i.e. for half-integer spin chain. This approach has been recently appliedlijiia to spin ladders where there is also 
a parity effect which is given by the number of legs of the ladder. Indeed, spin ladders with even number of legs 
are generically gapped while odd-numbered ones are gapless. This effect can be explained in the NLcrM framework 
with a topological term which is 9 — 2-ks x ni where ni is the number of legs. It is also of great interest to consider 
generalized spin ladders with various types of bond alternation, and/or additional exchange couplings. For example, 
a spin ladder with a diagonal coupling interpolates smoothlytj between the s=l chain and the two- leg spin ladder. 
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Previous investigationsE-3~t3 of these more general situations have used the NLerM but with the impossibility to treat 
cases including ferro bonds as in the ferro-antiferro chain introduced by Hida. 

In this work, we present a generalized NLerM that includes additional massive modes and we show that this model 
is able to reproduce the s=l limit of alternating chain of Hida. We discuss the construction of this new effective theory 
in the path-integral framework although it could also be done in the hamiltonian formalism with identical results. 

In Sect, [n] we treat Hida's chain case. We compute the spin-wave spectrum, motivate the NLerM approach and 
compute an approximate spin gap formula for the whole range of 7. In sect. [II, we apply the same formalism to 
White's mapping from the spin-1 chain to the antiferromagnetic two-leg ladder. Sect. IV contains our conclusions. 



II. FERRO- ANTIFERROMAGNETIC SPIN CHAIN 



A. Spin- wave spectrum 



A NLerM for the spin chain is a non perturbative field theory which is built on the lowest energy modes of the 
classical theory of the chain. We first study this alternating spin chain at the classical level. The equations of motion 
for the ferro-antiferromagnetic spin chain follow from -£rij = i[H, nj and are given by : 



^n 2l = -n 2 , x (-7112,+! + n 2l -i), 
^n 2l+ i = n 2i+ i x (-71124 + n 2l+2 ), 



(2.1) 



where is the classical spin vector at site i. Then we linearize those equations around the Neel order configuration 
n 2 i = (— l)'su + a, and n 2 j + i = (— l)*su + /3$ where u is an arbitrary unitary vector and at and /3j are orthogonal 
planar vector fields. Note that u breaks the natural 0(3) symmetry of the hamiltonian. We represent both ati and Pi 
by complex fields and Xi- The linearized equations of motion then reads : 



^ n = i a (-i)»(- 7Xn + Xn _ x + (1 + 7 )f„), 

f t Xn = is(-l) n (-rt n + Cn+l + (1 + 7)Xn). 

A plane wave Ansatz for those fields is : 

£„ = e^+ kn \a{k) + (-l) n b(k)), 
Xn = e^ +kn \c{k) + {-l) n d{k)). 

The eigenstates (a, b, c, d)(k) of energy u> and momentum k are the eigenvectors of the following matrix 



(2.2) 



(2.3) 





1+7 




-7 + e lk I + 7 

with eigenvalues ui{k). The energy w is a solution of its characteristic polynomial 



1 + 7 —7 — e 

- 7 + e- lfc 
-7 - e lk I + 7 

1 + 7 



— ik 



(2.4) 



(7)' 



47(1 + 7) [— J + 47 2 sin 2 k 



0. 



(2.5) 



The dispersion relation is given by : 



= 5^27(1 + 7) ±2 7v /(l+ 7 ) 5 



sin k. 



(2.6) 



The two lowest positive excitations correspond to momenta k c = and k c = it. Linearizing the dispersion relation 
around k c , we obtain for those modes the linear relation ui = s^J \d^\k — k c \. Hence the velocity of those spin waves is 

v = Syjjj^. They stand for the two Goldstone modes of the (classically) broken 0(3) symmetry. We will ground our 

NLerM on them, assuming they are slowing varying modes, in an theory where the rotational symmetry is unbroken. 
We also reasonably assume that the massive modes do not interfere significantly with the Goldstone modes, their 
energy being of the order of 2sy/j(l + 7). 
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B. The NLcrM mapping 



Here we briefly recall the derivation of the NLcrM for an antifcrromagnctic spin-s chain with s even. So consider 
a spin-s antiferromagnetic spin chain with s even, characterized by the hamiltonian TL = Jaf^i-^i+i- This spin 
chain can be mapped onto a NLcrM . The mapping can be performed within the lagrangian formalism if one uses 
coherent states representation for the spin operators. The discrete action is made of the exchange interaction terms 
and of the Berry phases W[ni] of the spin vector field 

S = JdtY^Js 2 ni. n i+ i +Y^s.W[ni\. (2.7) 



The coherent state vectors are expanded according to 



n 2 i = hi - s$2i 

n2i+l = hi+l + S$2i+lj 



(2- 



where the vector field <J>i (local staggered magnetization) is unitary and the field L satisfy nj.lj = 0. Those fields can 
be hinted at from the study of modes in the semi-classical spin chain. Now suppose L is of order of magnitude a, the 
lattice spacing. This representation of the spin chain allows for a continuous form for the action S, which depends on 
the fields l(x, t) and $(x, t). Because we supposed 1 is of magnitude a, the action is quadratic in 1 so that the field can 
be integrated out. Now recall that the spin magnitude s is even. Then the Berry phases of the spin coherent states 
partially gather to form a topological term 9 — 2irs, the effect of which is null in the action. The final effective action 
in the field <& is : 

S = J dtdx± Q($*) 2 - c(d x ^ , (2.9) 
where g is the coupling constant of the NLcrM : g — 2/s and c is the velocity c = 2Jafs. 



C. Effective NLcrM 



Let us consider now the alternating spin chain. The main difficulty we encounter in trying to represent the Hida 
chain is the fact that it is inhomogeneous, so that writing an adequate continuous action out of the discrete hamiltonian 
is not a trivial task. Indeed, the fundamental microscopic structure is a block of two spins (call them 1 and 2, say). 
It will turn out that it is better to choose pairs of spins naturally coupled by the ferromagnetic link, preferably to the 
antiferromagnetic one. This choice corresponds anyway to Hida's idea of pairing those spins s=l/2 to make them an 
effective spin-1 in the limit where 7 goes to infinity (large ferromagnetic coupling). 

The most rigorous way we could imagine to handle the continuous limit would be to introduce two coherent states 
vector fields, one for each of the two sites. Unfortunately, this would yield an intricate action, due to the appearance 
of several equally contributing massive modes. In particular, the Berry phases contribution of those coherent states 
would not be any more easily recognizable as a topological invariant. 

Now consider a two-block structure, made of the four spins S^, S^, S^+i, and S^+i- The pairs S^S^+i and 
Slo S^i+i are appropriate candidates to generate two NLcrM models. A first (incorrect) idea, sustained by our wish 
to go to the continuous limit, is to assume that they form the same NLcrM . This assumption would be correct in the 
limiting case 7 goes to infinity, but far too crude in any other case. 

To cure this, at least partially, we add one extra field Aj to the two semi-classical NLcrM slow modes L and It 
represents small quantum fluctuations, remnants of massive modes, that may bring about effective corrections to the 
NLcrM action. 

Accordingly, the coherent states fields are decomposed as : 

' = hi - s$ 2l + aA 2l 

n 2»+i = Wi + s9 2i +i - aA 2t+ i 
- n ii+i = hi+i + s$2i+i + aA 2i+ i. 

Here s is the spin magnitude and a is the lattice spacing. Note that this lattice spacing is of the length of a two-spin 
block. The amplitude of the quantum fluctuation is of the order of this lattice spacing. This assumption make the 
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problem a tractable one, since we do not need to consider derivatives of the field Aj when expanding the action. We 
enforce it by setting a as a prefactor of the fluctuating field. Note that the standard momentum field 1 is implicitly 
assumed to be of order a. ._. 

We intentionally chose only one fluctuation field, contrary to Senechal's schemed where the coherent state fields 
are decomposed on as many possible independent fluctuation fields. We now justify the choice of this particular field 
by means of a path integral reasoning. 

We suspect that some highly fluctuating paths are contributing in the action. Indeed, exchange couplings vary 
on the microscopic scale by a macroscopic amount, and do not behave smoothly with respect to the position. As 
a consequence some irregular paths might be energetically favorable. Yet a straightforward continuous limit of the 
action would not retain them since they are not spatially regular. That is why we should enforce some possible (non 
derivable) contributing fluctuation in the paths. Now let us see why we chose this particular field Aj. The fields L 
and <pi parametrize the variation of the path between the two-block pattern (2i,2i + 1), spatially indexed by i. So 
we need one field to represent the variation inside the two-spin blocks. In the block 2i, the coherent state vectors 
and n 2 ^ differ by an amount of 2aA^. Since the microscopic pattern is a two-spin block, we have then no choice 
than to make the coherent state vectors T^i+i an d n 2i+i differ by an amount of — 2aAi, because the variation on the 
scale of the two-spin block are already taken into account within the NLcrM fields. And this exhausts contributing 
infinitesimal fluctuations of the path. 

The action of the spin chain we get through the coherent states representation is of the form 

S = J dt^2s 2 n 2l . n 2l+ i - 7^ s 2 n 2 ,: + i.n 24+ 2 + ^ s.W^n,]. (2.11) 

i i i 

One can expand the coupling terms in the action, then goes to the continuous limit, which yields 

{ dr 

S c = dt— [-4a 2 (l + 7 )A 2 - 4Z 2 + 4a 2 sd x $.A - aV(9 x $) 2 ] . (2.12) 

We emphasize the fact that the ferromagnetic exchange terms are of course expanded from its aligned configuration 
contrary to the antiferromagnetic exchange terms which are expanded froixL.thc: Neel order. That is the reason why 
NLer models derived for antiferromagnetic ladder or alternating spin chainEjllj cannot be straightforwardly applied 
to the present cases : they are not built upon the same semiclassical configurations. Note that since the measure 
element dx is a two-spin block, it is equal to a. The Berry phases of the spins are of the form J dtdxSn.n A dtn, where 
5n is the spatial variation of the field n, and give 

r dx 

S b = dt—s(4l + 2asd x <f>).$ Ad t <!>. (2.13) 

Since we sum up the contributions for a double two-spin block, we must divide the whole sum by a factor 2. We 
end up with 

S = J dxdt -2(1 + 7)A 2 - 21 2 - 2sA.d x <P - ^s 2 {d x <p) 2 + J dxdt [s 2 (d x Q) A 9t$ + 2sl.<I> A , (2.14) 

where we have set a = 1 for commodity. We recover a spin-2s topological (9-term with 9 — 4ars. This has the 
consequence that the topological term does not contribute, so that the spin chain is likely to be gapped. 
We next integrate over the fluctuation fields, that is to say 1 and A. We then obtain the NLcrM action 



S = / dtdx 



2 \ t 1 2 l + x 1 
The standard parameters of this NLcrM are the coupling constant g and the velocity c given by 



(2.15) 



Hv^r and c= Vt?7 (2 - 16) 

Note that c perfectly matches the velocity v we found for the classical spin waves. When we make 7 goes to infinity the 
coupling constant g goes to 1/s and the velocity goes to s. These are the expected parameters for the 2s-NLerM . The 
Hida's chain is build up of spin one half so that for s = 1/2 g goes to 2 and c to 1/2 in units of the antiferromagnetic 
exchange coupling. Those are the parameters of a spin-1 chain of exchange Jaf/^- This is the limit found by Hida 
for the alternating chain. Let us recall why the magnitude of the exchange coupling should be so. Indeed in this 
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limit the ferromagnetic pairs are in a triplet state. Because of rotation invariance, the limiting hamiltonian can be 
written as an effective spin-1 chain. The coupling can only be quadratic S^.Sj or quartic (Si.Sj) 2 . The quartic term 
is excluded in this limit because it corresponds to second order perturbation theory in 7 _1 . The prefactor in front of 
the spin-1 coupling can be determined with the help of the calculation of a single matrix element, which is easily done 
on the all-spin-up configuration and yields J cff = Ja_f/4. So that in this limit, our NLcrM is consistent with Hida's 
argument. 



D. Energy gap evaluation in the large-TV limit 

To evaluate the spin gap, we compute the mass generated by the NLcrM in the limiting case of a large number of 
components for the field $(/V — > oo). In that limit a closed expression can be obtained for it, thanks to the large N 
saddle-point approximation. So we made $ a N-component field and rescaled it by a factor VTV". We then enforce the 
unitary constraint on $ by means of a conjugate field A, so that the unconstrained partition function is 



-S+i j dxdt\(<i> 2 -N) ^ 27) 



D<f>D\e 

The saddle point equation for the complete action can now be safely derived, and we may look for a constant 
solution for A that we will call im 2 . In the limit of large N, m 2 is the mass generated by the NLcrM . 

dk du> 1 

g gc 

In the limit 7 goes to 0, the chain is totally dimerized and the velocity c goes to 0. Because of this unusual feature, 
and also because we wish to derive results valid for a large range of 7, we will not resort to a radial cut-off in the 
euclidean space-time of the NLcrM . Rather we will first integrate on frequencies, then integrate on the momenta, with 
a large-momenta cut-off A. So that instead of the usual radial cut-off, we integrate over a strip in the plane (k, w) 
along the w-axis. The reason for this is that when 7 goes to 0, the prefactor eg -1 of k 2 goes also to zero whereas the 
prefactor (cc/) -1 of u> 2 remains constant so that large frequencies are more and more relevant and must not be cut off. 
In the process (which actually corresponds to the decoupling of dimers), we lost the space dimension of space-time. 
The first integration is over u> and gives 



A dk 1 



-A 47r ^Jk 2 + f m 2 
Then integrating over k we can extract the mass of the NLcrM 

A 



(2.19) 



9 sinh(27r S . M-) 



(2.20) 

In the standard derivation of the mass generated by the NLcrM (see for exa one would exchange the function 

hyperbolic sine for the function exponential. Indeed, in order for this computation to make sense we must have 
y/gc~ 1 m << A. Since in the NLcrM related to the spin-2s chain, y~gc~^ is finite, the (then meaningless) hyperbolic 
sine function can be replaced with the exponential function. Yet for our alternating chain, it can't be done since in 
the limit 7 goes to the argument of the function vanishes. Hence to encompass the full range of 7 we must retain 
the hyperbolic sine function. 

Now we can evaluate the energy gap : 



A s = ^gt.m = A * j=. (2.21) 

sinh^Trsy^) 

Whatever the spin magnitude s, the gap of an antiferromagnetic pair of spins is equal to the gap between the triplet 
state and the singlet, i.e. Jaf- Hence when 7 — > 0, we can determine that A goes to A/(2n). This statement allows 
us to determine the cut-off which appears to be 2tt (actually 2ir/a but we set a = 1). We can then write : 
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(2.22) 



We have obtained a self-contained, estimate of the chain-gap for the spin magnitude s . For the spin-1 chain we 
obtain Ai = si ^ - , which is 1.08 far from the numerically^ known 0.41. We expect a better result for the spin-2 
chain, closer to the "large spin limit". We obtain A2 = sin 8 1 ^ r 27r , which is 0.094 fairly close to the numericalljO known 
value of 0.085. As for the correlation lengths, we obtain £1 ~ 2 to compare with the numerically known £1^6 whereas 
we obtained £2 ~ 43 to compare with the numerical value £2 ~ 49 . 

The figures labeled 3 and 4 are drawings of the curves of the spin gap w.r.t. —7 in the range 7 S [0, 1] and — 1/7 ki 
the range 76 [1, 00 [. We did so in order to compare our results to Hida's presentation of his numerical computation!! 
Our result agrees qualitatively on the whole range of 7 including the limiting case 7 goes to zero. 

To check that our NLerM approach is still valid in the limit 7 goes to zero, we can also compute the correlation 
length. It can be read on the saddle point equation : £ = y^fm- Then we get : 

£ = — sinh ( 2irs, j — 
5 2tt V V 1 

We check that the correlation length goes to zero like s^/7 when 7^0, which is expected since at this point the 
chain is made of decoupled dimers. 



(2.23) 



III. AN ALTERNATING SPIN LADDER 

InRef. (0), S. R. White introduced several mappings that interpolate smoothly between spin ladders and a spin-1 
chain within the Haldane'spin gap phase. Here we treat one of these mappings. It consists in an antiferromagnetic 
spin ladder with an additional diagonal bond in every plaquette formed by legs and rungs (see figure 2). The exchange 
coupling associated to this bond is ferromagnetic (we will denote it as D) and does not introduce any frustration in 
the ladder. The hamiltonian is given by 

h = j J2 s:s; 1+1 i^sX + ^sX +1 (3.1) 

n;a— 1,2 n n 

where all exchange couplings are chosen positive. When D goes to we recover the usual antiferromagnetic spin 
ladder. Whereas in the limit —D goes to infinity, the pairs of ferromagnetically bounded spin are in a triplet 
state. The effective hamiltonian can then be expanded by rotational invariance in terms of spin-1 couplings. The 
effective antiferromagnetic coupling constant is then evaluated on any matrix element of the hamiltonian and gives 
J rff = (2J+ K)/A. So this limit corresponds to an antiferromagnetic spin-1 chain. Since in the process, we remain in 
the Haldane gapped phase, we may apply our scheme to obtain an estimate of the gap as a function of the diagonal 
coupling D. 



A. Effective NLerM 



The microscopic pattern of the ladder is composed of the four spins S^, S^, S^+i, and Sj i+1 forming a square and, 
apart from the diagonal bond, three bonds: one on each leg, and the last one on one of the two rungs closing the 
square (see figure 2). As for the Hida's chain case calculations are done on a doubled cell. 

Possible candidates to spin pairs forming a NLerM are the nearest neighbours on the same leg. Pairs of site linked 
by a rung contribute to the same NLerM . Then, we will need only one extra fluctuation field to describe of fluctuating 
paths inside the elementary cell, before going to the continuous limit. 

Accordingly the coherent states fields are decomposed as : 

' = hi - s$2i + aA 2i 
nf< = hi - s$2i + aA 2i 
n 2i+i = hi+i + s$2i+i 

> n 2*+l = hi+1 + S$2i+1 



aA 2 i+i 
aA 2i+ i 



(3.2) 



G 



In the following we set J = 1, K — p and D = — 5 so that the coupling constants are expressed in units of the 
longitudinal antiferromagnetic coupling J. One can then expand the coupling terms in the action, which yields : 

dt— \-4a 2 (p + 5)A 2 - 4(2 + p)\ 2 + Aa 2 5sd x 1>.A - a 2 (p + 5)s 2 (d x $) 2 } (3.3) 
a 

with the same care for the ferromagnetic bonds as was done previously. Like for the Hida's chain the measure element 
dx is a two-spin block and is equal to a. The Berry phases of the spins are of the form J dtdxSn.n A <9 t n and give 



dx 

dt— 2sl.$ A 

a 



(3.4) 



Since we sum up the contributions for a double two-spin block, we must divide the whole sum by 2. We end up 
with : 



S 



dxdt 



-2(2 + p)\ 2 - 2{p + 5)A 2 - 2sA.d x $ - ^s 2 (2 + 5)(d x $) 2 



dxdt [21s<& A 9t$] 



(3.5) 



where we have set a = 1. We next integrate over the fluctuation fields, that is to 1 and A. So that we finally obtain 
the NLcrM action 



S 



dtdx 



2 + p y ' 2 K 6 + p 



(3.6) 



The standard parameter of this NLcrM are 



• 9 = ;i 



2 + ^ 



and 



sJ(2 + p)(2 + 



5p_ 
S + p 



Now in order to stick to White's notations, we set p = 1 to get : 



1 3(1 + 6) 



s 3(2 + 35) 



s v 2 + 35 2 
These are the coupling and the velocity of the effective NLcrM 



1 



(3.7) 



(3.8) 



B. Energy gap evaluation in the large A^-limit 



A line of reasoning similar to Hida's chain treatment can be applied to this NLcrM . In the large N-component limit, 
we can derive a saddle point equation and obtain the generated mass. The energy gap is then given by 



a « . r 271". 

= A.c. exp( ). 
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(3.9) 



Note that contrary to the Hida's chain case, it is not meaningful to stick to the hyperbolic sine function, since it is 
not more precise than the function exponential. With the previously computed coupling constant g and velocity c, 
we obtain 



. s 3(2 + 36) 

A -2V^r exp 



-2ns\ 



2 + 35 



, 3 ( 1 + 5 )/ ' 

Since when 5 — > oo, we should recover 3/4 of the gap Ac of the spin-1 chain, we can rewrite it as : 



I 2 + 36 
3(1 + 5) 



exp 



2tts 1 



35 



3(1 + 5) 



(3.10) 



(3.11) 



Hence we can relate the spin gap of the antiferromagnetic spin-1 chain to the spin gap of the antiferromagnetic 
two-leg spin 1/2 ladder by the formula : 



Al 
Ac 



■ exp 



1 



(3.12) 



With the value Ac — 0.41, we obtain the estimate A/, ~ 0.45, close to the knowrc 2 ] value 0.50. In figure 5, we have 
drawn the curver©f the normalized spin gap A^/(3/4.Ac) w.r.t. 1 — 2/n. arctan(5. We can compare the curve with 
data from WhiteS-3. Not only does our result agree qualitatively, but it is also quantitatively quite good. 
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IV. CONCLUSION 



We have constructed non-linear sigma models appropriate to the description of the properties of some generalized 
spin chains including ferromagnetic exchanges. In the case of the alternating fcrro-antifcrromagnetic spin chain, our 
treatment reproduces correctly the smooth crossover from decoupled dimers to the Haldane phase. This approach 
may be of relevance to the study of the compound CuNt^Oe which is such an alternating chain and has a spin gapE3. 

We have also treated a ladder including ferro bonds so that the Haldane phase can be reached. Here again there is 
good agreement with numerical data. 
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FIG. 1. The alternating ferro-antiferromagnetic spin chain. The coupling 7 is the strength of the ferromagnetic bonds. The 
unit cell needed to construct a NLcrM contains four spins. 
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FIG. 2. A spin ladder with diagonal couplings S. The unit block to construct a NLcrM contains again four spins. 
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FIG. 3. The gap of the alternating chain as a function of 7. The point with decoupled dimers is on the right. 
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